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Abstract 

We study the charmful three-body baryonic B decays with D^*'^ or J/'f in the final state. We 
explain the measured rates of B'^ npD*'^ , B^ ppD^*^^, and B~ — > ApJ/^ and predict the 
branching fractions of B^ ApZ)*+, B^ T.^pD*+, B~ ApD^, and B'^ ApD*^ to be of order 
(1.2, 1.1, 1.1, 3.9) X 10~^, respectively. They are readily accessible to the B factories. 

PACS numbers: 13.25.Hw, 14.20.Pt, 14.40.Lb 
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I. INTRODUCTION 



There are several salient features of the charmless three-body baryonic B decays B BB'M 
with B being a baryon. First, a peak near the threshold area of the dibaryon invariant mass 
spectrum has been observed in many baryonic B decays Q, [s, li, [s], 6, Q, [s], [o], Q. The so- 
called threshold effect indicates that the B meson is preferred to decay into a baryon-antibaryon 
pair with low invariant mass accompanied by a fast recoil meson. Second, none of the two- 
body charmless baryonic B decays has been observed so far and the present limit on their 



branching ratios has been pushed to the level of 10 
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12] . This means that three-body final 



states usually have rates larger than their two-body counterparts; that is, r{B BB'M) > 
r{B BB'). This phenomenon can be understood in terms of the threshold effect, namely, the 
invariant mass of the dibaryon is preferred to be close to the threshold (for a review, see (l3j). 
The configuration of the two-body decay B — ^ BB' is not favorable since its invariant mass is 
ttlb- In B ^ BB'M decays, the effective mass of the baryon pair is reduced as the emitted 
meson can carry away a large amount of the energies released in baryonic B decays. Third, 
in the dibaryon rest frame the outgoing meson tends to have a correlation with the baryon or 
the antibaryon. Experimentally, the correlation can be studied by measuring the Dalitz plot 
asymmetry or the angular distribution of the baryon or the antibaryon in the dibaryon rest 
frame. For example, the proton in the dibaryon rest frame of B^ ppK~ is found to prefer 



to move collinearly with its meson partner [3| 
move back to back most of the time ] . 



whereas the proton and vr in i? 



While various theoretical ideas 113. 
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16|, 
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I20I have been put forward to 



explain the low mass threshold enhancement, this effect can be understood in terms of a simple 
short-distance picture 2l| . One energetic qq pair must be emitted back to back by a hard gluon 



in order to produce a baryon and an antibaryon in the two-body decay. This hard gluon is 
highly off mass shell and hence the two-body decay amplitude is suppressed by order of a^/g^. 
In the three-body baryonic B decays, a possible configuration is that the BB' pair is emitted 
collinearly against the meson. The quark and antiquark pair emitted from a gluon is moving 
nearly in the same direction. Since this gluon is close to its mass shell, the corresponding 
configuration is not subject to the short-distance suppression. This implies that the dibaryon 
pair tends to have a small invariant mass. 
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While the above-mentioned short- distance picture predicts a correct angular distribution 



pattern for the decays B ppn , B A^j»7r"'~, and B —>■ Apj [4, |8|, 122] , it fails to 



explain the angular correlation observed in B^ ppK~ js], [sl and B~ — ^ Apvr" The 
intuitive argument that the K~ in the pp rest frame is expected to emerge parallel to p is not 
borne out by experiment Likewise, the naive argument that the pion has no preference 

for its correlation with the A or the p in the decay B^ Apn^ is ruled out by the new Belle 
experiment j?! in which a strong correlation between the A and the pion is seen. Therefore, 
although the short-distance picture appears to describe the global features of the three-body 
baryonic B decays, it cannot explain the angular correlation enigma encountered in the penguin- 
dominated processes B^ ppK~ and B~ Apir^ . 

The study of the charmful baryonic B decays B BB'Mc may help improve our under- 
standing of the underlying mechanism for the threshold enhancement and the angular distri- 
bution in three-body decays. First, the aforementioned three features also manifest themselves 



in B ^ BB'M^. An en 

the decay B ppD^*^ 



lancement at the low dibaryon mass has been seen, for example, in 
2J], while the Dalitz plot oi B ppD^*^ 2J] with asymmetric 
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distributions signals a nonzero angular distribution asymmetry. Therefore, B — * BB'Mc is a 
mirror reflecting the phenomena observed in B ^ BB'M. Second, since the decay rates of 
B — * BB'Mc are 10-1000 times larger than that oi B BB'M, in principle this will render 
the detection of angular distributions and Dalitz asymmetries in the latter easier. Theoreti- 
cally, the relevant topologic quark diagrams for charmful decays are simpler than the charmless 
ones as the latter receive not only tree but also penguin contributions. In this work, we shall 
therefore focus on B BB'Mc with Mc = D^*^ or J/\l/, and our first task is to see if we can 
explain the branching fractions. 

This paper is organized as follows. The experimental data for the decay rates of i? ^ BB'Mc 
are collected in Sec. II. The formulism is set in Sec. Ill with special attention to various 
transition form factors. We then proceed to have a numerical analysis in Sec. IV and discuss 
some physical results in Sec. V. Section VI contains our conclusions. 
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II. EXPERIMENTAL DATA 



In this section, we summarize the measured branching ratios for the three-body charmful 
baryonic B decays B BB'Mc: 



Br{B^ npD* 



BriB 







Br{B^ ppD 



*0\ 



Br{B^ 
Br{B- 
BriB- - 



>ApD; 

^ ppD' 
ppD*' 



= (U.btli ± 2.7) X 10"^ [CLEO] , 
= (1.13 ± 0.06 ± 0.08) X 10-^ [BaBar], 
= (1.18 ± 0.15 ± 0.16) X 10^^ [Belle], 
= (1.01 ± 0.10 ± 0.09) X 10^^ [BaBar], 
= (1.20j:[!:|±0.21) X 10"^ [Belle], 
= (2.9 ± 0.7 ± 0.5 ± 0.4) X 10"^ [Belle] 

< 0.15 X 10"^ (90% C.L.) [Belle] , 

< 0.15 X 10"^ (90% C.L.) [Belle] , 



for Mc = y, y, y, bej, and 



Br{B' 



ApJ/^) = (11.6^^:^1^:8) X 10^^ [BaBar] 



;ii.6± 2.8:^2:3) X 10" 



[Belle] , 

Br{B- TPpJ/-^) < 1.1 X 10"^ (90% C.L.) [Belle] , 
Br{B^ ^ ppJ/<U) < 1.9 X 10"^ (90% C.L.) [BaBar] 
< 8.3 X 10"^ (90% C.L.) [Belle] , 
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for = J/ijj 
by CLEO in 2001 
in rates. 



(1) 



(2) 



Among various measured decay modes, B^ npD*^ was first observed 
25| . Note that the decays B^ ppD^ and B^ ppD*^ have similar results 



III. FORMALISM 



The effective Hamiltonian responsible for charmful baryonic B decays reads 29] 

Heff = ^V,tV;^,{ciOi + C2O2) + H.c. , 
v2 



(3) 
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FIG. 1: Two types of the B BB'Mc decay process: (a) current type and (b) transition type. 

where Vij are the CKM matrix elements and q the Wilson coefficients. The four-quark operators 
Oi are defined by 

Oi = iq'^qa)v-Aici3bf3)v-A , O2 = {qMv~Aicf3ba)v-A , (4) 

where {q'^qj)v-A denotes ^■7^(1 — 75)'?j with q' = d, s, q = u, c, and i, j the color indices. 

To proceed, we shall adopt the generalized factorization approach 30, iMl in which the vertex 
and penguin corrections to hadronic matrix elements of four-quark operators are absorbed in 
the effective Wilson coefficients c^^ so that the fi dependence in transition matrix elements is 
smeared out: 



J]Q(/x)(Q.(/i))^5^cf(g.)viA, 



(5) 



where the subscript VIA means that the hadronic matrix element is evaluated under the vacuum 



insertion approximation. This a ppr oac. 



1 has been wel 



baryonic B decays 
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19 



32 



33 



34 



35 



36 



37 



applied to the stud y o f three-body 
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39 



40 



42 



44| . 



Under the factorization approximation, the decay amplitudes can be classified into three 
different categories: the current-type (class-I), the transition-type (class-II), and the hybrid- 
type (class-Ill) amplitudes. The class-I current-type amplitudes which proceed via a color- 
allowed, external I^-emission diagram as depicted in Fig. [T](a) are given by 



^^=K.KV«r'(BB'|(g'n)v._A|0)(/^(*)|(c6)v^_A|i?) 



V2 



(6) 



with af'*' to be specified later. By fixing q = u to avoid the charm baryon production, all 
possible decay modes in this category are 

5° jnp , S^A , AS+ , S^SO , , S^So} D^*^+ for q' = d , 

'Ap , S°p , S-n , S-SO , S-A , forg = s. (7) 



5^ 



The class-II transition-type amplitude via the color-suppressed internal W emission diagram 
[Fig. [11(b)] reads 

At{B^BB'M,) = ^KbV^Vaf (M,|(cg)v_A|0)(BB'|(g'6)y-A|i?), (8) 

with 02^'' to be given later. The possible decay modes in this class are 

jnp , S"A , AS+ , S^SO , , E'^^ J/^ for q' = d , 

B- , TPp , S-fi , E^TP , , J/^ for g' = s , (9) 

for the charged B system, and 

fpp, AA, S+S+ , S-S^ , ASO 1 , 
1^ , S^A J 

B° |An , , S+p , S^SO , E^A , S-S^ j J/^) for q' = s , (10) 

for the neutral B system. The class-Ill hybrid-type amplitudes which consist of both the 
current-type and the transition-type transitions are given by 

AniB ^ BB'Z}(*)) = ^K.Kv{«r'(BB'|(g'«)y_^|0)(DM|(c6)v^_^|5) 

+af *'(DW|(c«)y_^|0)(BB'|(g'6)^_^|5)| . (11) 

The allowed modes are 

B- 

B^ 



jnp , S-A , AS+ , S^SO , , for g' = , 

|Ap,SV,S^fi,S~SO ,H^A,S°E+|dW° forg' = s. (12) 



We now can have a qualitative understanding of the data in ([T]) and ([2]). The nonobservation 
of B~ ppD^*^^ is due to the fact that it proceeds via h — > ucd at the quark level. This leads 
to a suppression of \VubV*^/VcbV*^\'^ ~ 10^"^ compared to its neutral partner. Likewise, it is 
expected that Br{B^ ppj/"^) = \Vcd/Vcs\'^ Br{B' ApJ/^) ~ 10"^, consistent with 
the experimental upper bound for this decay mode. As for the comparable decay rates for 
B^ ppD*^ and B^ ppD'^, it has to do with the baryonic form factors, which we are going 
to elaborate on later. 
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A. Decay constants and form factors 



To proceed, we need the information of the decay constants and form factors relevant for 
charmful baryonic B decays. Decay constants are defined by 

(K|gi7M92|0) = mvjv.el , (13) 

for the pseudoscalar and the vector meson Vc- The B to D^*'^ transition form factors are 
parametrized as [45 1 



(D|c7^6|5) 



(Pb+PdT- 



m 



B 



m 



D 11 



m 



B 



m 



rriB + 



(D*|C7;.75&|5) = I 



■ q 



— I 



{PB+PD*)^^ - 



m 



B 



m 



D* 



is* ■ q) 



it) 



TTlB + mo* 



(14) 



where t = with q = Pb — Pd* = Pb + J'g'- 
For the dibaryon creation, we write 

(BB'|gi7^g2|0) = u{p^){ F^{t)^ ^ + 



me + mi 



u(pb) [Fi(t) + F2(t)]7;.+ 



me + mi 



-iPB'-PB)^ >vip 



(BB' 1 917^75^2 1 0) = u{PB)<9Aithfi + 



hAit) 



9m \l5vipB'] 



(15) 



m-B + mg, 

where u{v) is the (anti-)baryon spinor, and F12, Qai are timehke baryonic form factors. 
Note that there are two additional form factors in the form of ug^f and ua ^j^qp'-^^v . However, 
since we assume SU(3) flavor symmetry, we can neglect these two form factors as they vanish 

for conserved currents. The asymptotic behavior of form factors is governed by the pQCD 

counting rules {46, 



48|. In the large t limit, the momentum dependence of the form factors 



47 



Fi{t) and QAit) behaves as 1/t^ as there are two hard gluon exchanges between the valence 
quarks. More precisely, in the t ^ 00 limit 



Flit) 



C 



Fi 



In 



A2 
^^0 



9Ait) 



a 



9A 



t2 



In 



A2 
^^0 



1 -7 



(16) 



TABLE I: The parameters Cp^ and C^^ in Eq. ()16p . expressed in terms of the combination of C\\ and 
Cm-, where the upper (lower) sign is for Cp^ (Cg^). 



^ 







Ci^i(Cg.4) 


np 


1^1 1 ^ 5^11 


Ap 


-^/i^ii 


S-A 
AS+ 


i^^ll ^ ^'^^ 


TPp 
S~n 

H-A 


3 ^'11 + a*-"!! 



where 7 = 2 + 4/(3/5) = 2.148 with /5 being the QCD p function and Aq = 0.3 GeV. In the 
asymptotic t ^ 00 hmit, both F2{t) and /i^(t) have an extra 1/t dependence relative to Fi 
and qa owing to a mass insertion at the quark line 38, Q, [s^. However, the form factor is 
related to qa by the relation 



hA 



(me + mg, 
t 



-9 a 



(17) 



through the equation of motion. Hence, in ensuing numerical analysis we will keep /lyi(t) and 
neglect F2{t). Under the SU(3) flavor and SU(2) spin symmetries 48|], the parameters Cp^ and 
Cg^ appearing in ^ BB' transitions are no longer independent but are related to each other 
through the two reduced parameters C\\ and C-^ (see Table H]). 

For the three-body transition B —>■ BB', its most general expression reads 

{BB'\q\b\B) = 

iuiPB)[gil^t + 92icr^,^p'' + g^p^ + g^{p^, + Pb)^ + 5-5 (Pb' " Pb) ^,]l5V {Pb') , 
(BB'|g'7^756|i?) = 

^^(Pb)[/i7m + f2i(T,,uP'' + hPi, + /4(Pb' + PB)^, + h{Pw - Pb)^v{p^,) . (18) 

with p = Pb — Pb — Pb'- principle, the form factors fi and Qi depend on not only the 
invariant mass t of the dibaryon, but also the invariant mass (i.e. the variable u or s) of 
one of the baryons and the emitted meson. Indeed, such a momentum dependence has been 
studied in the framework of the pole model in which some intermediate pole contributions to 



3, Q, Q, 



the three-body transition of -B ^ BB' are considered pj, iizi, m Il9l. |32|]. However, since the 
momentum dependence of the transition form factors on the variable s or u is poorly known, 
one approach which is often employed in the literature is that one first parametrizes them in 
the form of a power series of the dibaryon invariant mass t and assumes that the dependence 



3, 



38 



39 



4G 



41 



5l|. 



on other variables may be lumped into the coefficients of 1/t expansion 
This is the approach we will adapt in this work. Since three gluons are needed to induce the 
B BB' transition, two for producing the baryon pair and one for kicking the spectator quark 
in the B meson, pQCD counting rules imply that to the leading order 



9^it) 



'9i 



t3 



(19) 



Just as the previous case for vacuum to the dibaryon transition, under the SU(3) flavor 
and SU{2) spin symmetries, the parameters Dg-^, Df-^, Dg^ and Df^ can be expressed in terms 
of the reduced parameters D\\, D-^ and (see Table HTj) . and the derivation is presented in 
the Appendix. Interestingly, the decays B^ S"''S+ and S^H" are prohibited by SU(2) spin 
symmetry. 

Since many of the three-body baryonic B decays show an enhancement at the low dibaryon 
invariant mass squared t, low values of t produce the dominant contributions to these decays. 
Recently, BaBar has measured the phase-space corrected dibaryon invariant mass distributions 
for various B D^*^pp decays 2J]. This amounts to measuring the squared amplitude \A\'^ 
versus t. It turns out that the t dependence of \A\^ can be represented sufficiently by a single 
1/t" behavior. Therefore, it is justified to apply the pQCD counting rules valid at large t to 
form factors at small values of t. 



IV. NUMERICAL ANALYSIS 



We need to specify various input parameters for a numerical analysis. For the CKM matrix 
elements, we use the Wolfenstein parameters A = 0.807± 0.018, A = 0.2265 ±0.0008, p = 0.141 



and T] = 0.343 [5^. For the decay constants, we use (53|, [54 1 



Ud, fD',M = (0.22, 0.23, 0.41) GeV 



(20) 
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TABLE II: The parameters Dg^, Dj^, Dg-, and Dj. in — > BB' transition expressed in terms of the 
reduced parameters Dw, Df, with i = 2,3,4,5, where the upper (lower) sign is for Dg^ (-^/i)- 



B- 



np 

AS+ 



B 







pp 1 

nfi 1 

AA 

=0=0 J 



V6 



II V6 II 



^^11 

-sjl^l 



3V2 



Dtt 



1 



^11 ^ i^H 



5^11 T 1^^ 
l^ll T ii?^ 
I^ll =F il?^ 







3^11 



3 

2v^ 



^11 
^11 



B- 



Ap 



l^ll 



i^ii 



^^11 







An 



■373^1 
1 



3^2 



3^11 ± i^ii 



^^11 



3v^ 



^^11 ^ 75^F 



l^ll 



3V2'^II ^ 3v^^ll 
3-^11 =^ 3^11 

"373^11 ^ 373^1 
75^11 ^ 75^F 



-Di 

2 



373^fl 
Te^fl 
i^ll 



3 r)i 
2^11 



3^2^11 

-Di 

3 II 
3^2^11 
V6 II 

-Di 

3 



Following 55|, the momentum dependence of i? — Z^*^*^ transition form factors in Eq. is 
parametrized as 

m 



b(q^ /mf '^'^ ' 



We shall use the parameters a, h and -F(O) obtained in 55| ( see Table UTTll . 



(21) 



TABLE III: Various B — > D^*^ form factors taken from 



55|. 





Fi Fq 


Vi Ao Ai A2 


m 

a 
b 


0.67 0.67 
1.25 0.65 
0.39 0.00 


0.75 0.64 0.63 0.61 
1.29 1.30 0.65 1.14 
0.45 0.31 0.02 0.52 



For the parameters C\\ and Ctt in Eqs. ( IT5l [T6|) and Table [U we use the data of e+e 



56 



pp, nn 



571] to determine their magnitudes and the decay rate of — > npD*^ to fix their relative 
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sign 



(q,, C-rr) = (67.9 ± 1.4, -216.9 ± 23.5) GeV^ , 



(22) 



with x^/d-O-f = 1.6. As for the parameters Du and D-n in Eqs. f lTSj and Table [TTl we 



em 



qy the observed rates of ppD^, B ppK* , B^ ppK*^, and B 



mode to obtain 
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PPTT' 



24| in conjunction with the measured angular distribution of the last decay 



(L>ll, L>^) = (67.7 ± 16.3, -280.0 ± 35.9) GeV^ 
{Dl Dl Dl Dl = 

(-187.3 ± 26.6, -840.1 ± 132.1, -10.1 ± 10.8, -157.0 ± 27.1) GeV^ , (23) 



with x^/d-o.f = 0.8. To calculate the decay rates, we use the relation 

1 lAP 



dV 



-dm-g-Qidm-Q/j^j 



(24) 



(27r)3 32M| 

with m^g, = (Pb + ^b')^ ^B'M — (-^B' ~l~ -Pmc)'^: and the baryon spins have been summed 
over in The numerical results for the branching ratios are summarized in Tables HVllVIIt 
where the theoretical errors come from the uncertainties in form factors. 



V. DISCUSSION 



In the generalized factorization approach, the coefficients ai and a2 appearing in the decay 



amplitudes Ac, At and An are given by 29 



30 



3l| 



ai = cf + \ + cfxi 



02 



„efr 



„eff 



(25) 



with (c^ , d 



(1.169, —0.367). Since the nonfactorizable effects characterized by the 



d^^Xj terms are not calculable, we will fit them to the data. We see that 02 is dominated by the 
nonfactorizable term due to the large cancelation between and df^ /?>. Hence, nonfactorizable 
contributions to color-suppressed modes are sizable. Factorization works if the parameters ai 
and 02 are universal; namely, they are channel by channel independent. Empirically, this is 
supported by the experimental measurements 58|]. Two-body mesonic B decays suggest that 
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TABLE IV: Branching ratios for the current-type — > BB'D^*)^ decays. The results for D* pro- 
duction are shown in parentheses. 



BB' 


Br{B^ BB'L>(*)+) X 10^ 


BB' 


Br{B^ BB'D^*')+) x 10^ 


np 


5.2 ±0.5 


(14.5 ± 1.4 ) 


Ap 


3.4 ±0.2 


(11.9 ±0.5) 


S^A 


0.3 ±0.1 


(1.3 ±0.2 ) 




2.8 ±0.7 


(10.5 ±2.7) 


AS+ 


0.3 ±0.1 


(1.3 ±0.2 ) 




5.5 ± 1.4 


(20.9 ±5.4) 




0.05 ±0.01 


(0.23 ± 0.02) 




0.6 ±0.1 


( 3.0 ±0.3) 




0.05 ±0.01 


(0.23 ±0.02) 


E^A 


0.2 ±0.1 


( 1.3 ±0.4) 




0.07 ±0.02 


(0.4 ±0.1 ) 




1.2 ±0.1 


( 6.2 ±0.6) 



TABLE V: Same as Table HV] except for the transition-type S° BB'L>(*)° decays. 



BB' 


Br{B^ BB'L>(*)0) x 10^ 


BB' 


Br{B^ BB'L>(*)°) X 10^ 


pp 


1.1 ±0.2 


(1.0 ±0.4 ) 


An 


8.2 ±2.7 


( 9.8 ±3.7) 


nn 


3.9 ±1.3 


(5.8 ±2.3 ) 


S% 


0.7 ±0.1 


( 0.6 ±0.2) 


AA 


0.02 ±0.01 


(0.03 ± 0.02) 


S+p 


1.4 ±0.3 


( 1.2 ±0.4) 




0.10 ±0.03 


(0.07 ±0.03) 


HOEo 


1.2 ±0.4 


( 0.7 ±0.3) 




0.4 ±0.1 


(0.3 ±0.1 ) 


HOA 


0.09 ±0.04 ( 0.10 ±0.04) 




0.022 ±0.004 (0.012 ±0.005) 


e^tF 


2.2 ±0.7 


( 1.4 ±0.5) 




0.18 ±0.05 


(0.12 ±0.04) 









TABLE VI: Same as Table HV] except for the hybrid-type B' BB'Z)(*)° decays. 



BB' 


Br{B^ BB'L>(*)°) X 10^ 


BB' 


Br{B- BB'i:)(*)°) X 10^ 


np 


13.0 ±2.9 (46.6 ±5.3) 


Ap 


11.4 ±2.6 (32.3 ±3.2) 


S-A 


0.6 ±0.1 (2.6 ±0.3) 




3.3 ±0.6 (15.2 ±2.4) 


AS+ 


0.6 ±0.1 (2.6 ±0.3) 




7.0 ±1.2 (30.2 ±4.7) 




0.3 ±0.1 (0.6 ±0.1) 


E-TP 


1.9 ±0.3 ( 6.2 ±0.5) 




0.3 ±0.1 (0.6 ±0.1) 


E-A 


0.5 ±0.1 ( 1.2 ±0.2) 


e~eP 


0.09 ±0.01 (0.6 ±0.1) 




3.7 ±0.7 (12.7 ±1.1) 
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TABLE VII: Branching ratios for the transition- type B — > BB' J/^' decays. 



BB' 


Br{B~ BB'J/^) X 10^ 


BB' 


Br{B- 


BB' J/l-) X 10^ 


np 


14.4 ±3.2 


Ap 




11.6 ± 2.9 


S-A 


0.04 ±0.02 






0.13 ±0.02 


AS+ 


0.04 ±0.02 


S~n 




0.24 ±0.05 




0.29 ±0.09 






22.9 ± 5.8 




0.29 ±0.09 


E-A 




0.52 ±0.33 




0.37 ±0.08 






44.8 ± 10.8 



TABLE VIII: Branching ratios for the transition-type B^ — > BB' J/^* decays. 



BB' 


Br{B^ BB'J/^-) X 10^ 


BB' 


Br{B^ BB' J/^) X 10^ 


pp 


1.2 ±0.2 


An 


10.4 ±2.6 


nn 


6.7 ± 1.7 


YPn 


0.11 ±0.02 


AA 


0.0010 ± 0.0004 


S+p 


0.24 ±0.05 




0.13 ± 0.04 


E^TP 


20.6 ±4.8 




0.54 ± 0.16 


E^A 


0.5 ±0.3 




0.35 ± 0.08 


e-yF 


42.1 ± 10.4 


rAEO 
tsOA 


0.09 ± 0.02 







ai ~ C(l) and as ~ 0(0.2 - 0.3) |58|, l59fl. It will become plausible if af= fitted from the 
baryonic B decays lie in the aforementioned ranges. Indeed, the results given by 



a?* = 1.23 ±0.19 , a?* 



0.33 ±0.04, a2^* 



0.17 ±0.03 



(26) 



agree with the above expectation. Hence, we shall assume the validity of factorization in 
charmful baryonic B decays. 

There are several salient features we learn from Tables HVTl VII (i) B{B'^ BB'D*+) / B{B'^ — *• 
BB'i5+) > 3 for class-I (current-type) decays; (ii) for class-II (transition-type) decays, B{B^ — >• 
B'B'D*^)/B{B^ B'B'D^) > 1 or < 1, depending on the modes under consideration, and (iii) 
B{B^ BB'D*°) /B{B^ BB'D^) > 3 for hybrid-type decays except for BB' = E~TP, 
and S'A. 
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To get insight into the above-mentioned first feature, we take BB' = np as an example. In 
the dibaryon rest frame we have pb = —f>w ^"^^ Pb = Pdm- the branching fractions are 
dominated by the threshold effect, we can consider the dibaryon in the nonrelativistic limit 
where = (pe + Pb')i^ ~ {2m, 0, 0, 0) with E-b{E-q,) ^ m. Then we find 

(BB'|gi7^g2|0) ^ -2m (0, 1, ±z, 0)Fi(4m2) , 
(BB'|gi7^75g2|0) ^ -2m(l, 0, 0, 0)(^7^(4m2) + /i^(4m2)) . (27) 

However, the vacuum to BB' transition induced by the axial-vector current is suppressed as 
/i^(4m^) ^ — 5f^(4m^) followed from Eq. (fT7|) . From Eqs. (j6]) and (fT4|) we see that the decay 
amplitudes npD*^ and npD*^ are governed by the form factors Af^* and 

respectively. Thus the ratio of these amplitudes is given by 

A{B^ npD*+) _ me + mp* Ep* 

A{BO^npD+)^ 2mp \pp\' ^ ^ 

where we have assumed af* = af and taken the approximation of Af^*/Ff^ ~ 1 in- 
ferred from Table [XTTl Note that since numerically (m^ + mp*)"^ /{Am'jj) ~ 3, we thus have 
npD*^)/B{B^ npD^) ~ 3 as the ratio Ep*/\pp\ is close to 1. For heavy 
dibaryons, the ratio B{B^ BB'D*+)/i3(i?° — > BB'/)''') will become larger, for example, 
B{B^ E-^D*+)/B{B^ E-^D+) ~ [{niB + mp,)Ep,/{2mp\pp\)]'^ = 4.2. Among the 
current-type modes listed in Table llVt we notice that B^ ApD*'^ and B^ TPpD*^ have 
the largest rates 

B{B^ ApD*+) = 1.2 X 10"^ , 

B{B^ TPpD*^) = 1.1 X 10"^ . (29) 

A search of them at B factories is strongly encouraged. Note that the experimental studies of 
the angular distributions in the charmful cases may help solve the angular correlation puzzle 
in B^ Apn^ owing to the same BB' transition form factors appearing in both cases. 

Similar to the charmless cases of B^ ppK*~ and B^ ppK*^ 40|, |4ll], the decay rate of 
the class-II decay B^ ppD*^ is dominated by the form factors /2 and g2 defined in Eq. f[T5]l . 
Contrary to the D* production, the decay B^ — > ppD'^ receives contributions from the form 
factors /a and accompanied by m|,. Considering the dibaryon in the nonrelativistic limit. 
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gi\PD* \ + fiED*sm6p + 2g2mD*iV2\pD*\cos6p + mD*sm6p) 






qiEd + /i|pD|sin6'p + g-im\) 



the ratio of the amphtudes reads 

-I- fi Rr^ieinfl -I- 9 nr,m I ^/9 1 rTn* I pnefl -I- m einfl i 

(30) 

This ratio is expected to be of the order of unity since gi and /i terms in Eq. (!30l) are of the 
same order, and so are g2 and g^ terms. 

Class-Ill decays will have rates larger than class-I and class-II ones if the relative phase 
between Oi and 02 is small so that their interference is constructive. To proceed, we use af 
and extracted from npD*^ and ppD*'^, respectively, and neglect their phases. 

The results are summarized in Table IVI[ For example, 

B{B^ ApD^) = 1.1 X 10-^ 

B{B^ ApD*^) = 3.2 X 10"^ , (31) 

which are larger than S(-B° — > ApD~^) and B{B^ — *■ ApD*~^), respectively, by a factor of 3. It 
will be interesting to test this experimentally. 

As seen in Table Ivn] and IVIIIl for the baryonic decays with a J/ip in the final state, our 
prediction B{B~ TPpJ/'^) = 1.3 x 10^'^ is consistent with the Belle limit, 1.1 x 10"^, 
and B{B^ ppj /'^) = 1.2 x 10^^ is in accordance with the BaBar limit but slightly higher 
than the upper bound set by Belle [see Eq. ([2]) for the data]. We also find that the decays 



B^ (H^SO, S-A, S^SO, S^S-)J/* have branching fractions of order 10"^ 

As for the threshold peaking effect, while it manifests in the decay B^ ppD^*^^ the 
data clearly do not show the threshold behavior in B^ ApJ/"^ (see Fig. [2]). This can be 
understood as follows. In the latter decay, the invariant mass mj^p ranges from 2.05 to 2.18 GeV, 
which is very narrow compared to the nipp range in the 5° — ^ ppD^*^'^ decay. Consequently, the 
invariant mass distribution of dT/dm^p is governed by the shape of the phase space due to the 
relative fiat dependence within the small allowed m\p region. 

Finally we notice that since the transition-type decays share the same B BB' transition 
form factors with B^ PpT^^ , B~ ppK~ and B~ Ap7, measurements for distribu- 
tions/Daltiz plot asymmetries are useful for solving the puzzle of the angular distributions in 
the charmless baryonic B decays. For example, a recent Dalitz plot analysis of B^ ppD^*''^ 
by BaBar 2J] shows a possible indication that 

Ae{B° ppD^) > AeiB^ ppD*^) > , (32) 
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m 



2.5 3.0 
pp {GeV) 



2.5' 3.0 
mpp(GeV) 



2.00 



2.05 2.10 2.15 



2.20 



FIG. 2: Dibaryon invariant mass distributions for 
respectively. Experimental data are taken from 0,13. 



ppD^, ppD*^ and B 



ApJ/^, 



where the angular asymmetry Aq is defined by [5ll 

-+1 dr 



A, 



fn cos e - f\ -r^d cos e 

JO a cos J — 1 rf cos 



C^TT^d COS 9+ r° ;v^cicos^ ' 

JO a COS J — 1 a cos 

which is equal to (A^+ — ) / (A'^-i- + ) , where N± are the events with cos ^ > and cos ^ < 0, 
respectively. The above relation fl32p can be explained in our model: in the previous study 
5l| we have predicted that \Ag{B~ ppK~)\ > \Ag{B~ ^ ppK*^)\ and \Ag{B~ 



(33) 



m 



41 



pp7r~)| > \Ae{B~ PPP~)\- In our model, the charmful baryonic B decays have similar features 
of the angular distribution as the charmless cases. Therefore, we believe that our model together 
with the choice of the form factors will help understand the angular correlations observed in 
these decays. 



VI. CONCLUSIONS 



Within the framework of the generalized factorization approach, we have studied three dif- 
ferent types of charmful three-body baryonic B decays with D^*'' or J/^ in the final state. We 
found that the B BB' form factors extracted mainly from the data of charmless baryonic 
B decays can be used to explain the measured rates of B^ ppD^*^^ and B~ ApJ/\&. 
The branching fractions of B^ ApD*+, B^ TPpD*+ , B^ ApD^, and B' kpD*^ 
are predicted to be (1.2, 1.1, 1.1, 3.9) x 10~^, respectively. They are readily accessible to B 
factories. 
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Appendix: Relations for B BB' transition form factors 



Form factors of (B'B'\q'-f^b\B) and (BB'|g7^75&|i?) can be related to each other by the SU(3) 
flavor and SU(2) spin symmetries in the large t limit {t — ^ oo). As the currents V"^(A^) = 
57^(75)6 are used to create the dibaryon and annihilate the B meson with \B) ~ |&75'?'|0), we 
have 



{BB'\V^ -A^,\B) 



{BB'\2qLlMbLqR-bnq'M 

{BB'\2qL-fM + T^b)qR - '2qLlM + ^b)qL\0) 

(BB'I - 2mbqLlf.q'L + 2gL7/. MrI^) 

{bbVmIo) + (bbv;|o), 



(34) 



where = —2mbqLlfiq'L and = 2qL'y^ l/bq'R- Note that ( J^) is the source of the chiral-even 
(chiral-odd) current. In terms of the crossing symmetry, the antibaryon in the fina^ 



be transformed as the baryon in the initial state. Then we can follow Refs. 
parameterize the B' ^ B transition form factors as 

I-75 



37 



state can 



38 



48j to 



(B|J^|B') = 2zmfcU7^ 
(B| j;|B') = 2zu7^ 



75 

2 

1 + 75 



D- 



-D- 



u . 



75 



D'- 



u . 



2 2 

The above chiral-even and chiral-odd form factors and D'^ are given by 



(35) 



ef;Du + e±Du 



(36) 



where the amplitudes in which the currents interact with a valence quark with 

helicity parallel (antiparallel) to the helicity of the baryon. The parameters ejj" and ei are 
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defined by 



Y,{B,h\Q{{) + QU)\B',h') 

i<j 

J2{B,h\Q{k)\B',h'), 



(37) 



where the charge Q = Jq with i,j,k = 1, 2, 3 numbering the valence quarks in the baryon to be 

acted. Therefore, et,TTs is the sum of the electroweak charges carried by valence quarks in the 

11(11) 

baryon with hehcities parallel (antiparallel) to the baryon's helicity. Besides, the superscript 
+ (— ) denotes h' =| (|), the helicities of the initial baryon. 

Take 13^ pp as an example. With the helicity state of the proton given by 



b,T) 



\pA) 



we obtain 



b,nT) + b,Tn) + b,iTT) 

(|p,t) with t^i), 



i<j 

J2{p,UQi^) + Qi3)\pA) = o, 

i<j 

J2{pA \Qik)\pA) = , 

k 

Y,{pA\QmpA) = \. 



for the chiral-even J^, and 



I](p,ilQ(0 + Q(j)b,i) = o 



$^(p,TIQ(0 + Q(j)b,T) 

Y.^Pd\Qmpd) = ^. 

k 

$^(p,TIQ(A:)b,T) = o, 



(3^ 



(39) 



(40) 



for the chiral-odd J'^. When the sums of the electroweak charges ejy^||^ derived in Eqs. (139| l40l) 
are put back into Eq. (l35ll . (B| + J^|B') is given as 

+ ^M|m4^D|| + ^D^]7^ + [^D[|]7,,x/,|m. (41) 

This gives the form factor relations for (BB'|g7^6|i?) and (BB' 1^7^7561 5), which is given by 

12 12 
91 = 3^1, -3%,/. = 35,1 + 3%. 

<li = -n^-\D\, (i = 2,..-,5). (42) 

By the same token, we find the relation of gi = —fi holds in all B BB' transition form 
factors. 
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